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Abstract 

A picture is presented in which standard physics and its extensions are obtained from statistical 
counting and stochastic fluctuations, together with the geography of our particular universe in D 
dimensions. The inescapable predictions include super symmetry, SO(d) grand unification, Higgs- 
like bosons, vanishing of the usual cosmological constant, nonstandard behavior of scalar bosons, 
and Lorentz violation at extremely high energies. 
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I. INTRODUCTION 



For a theory to be viable, it must be mathematically consistent, its premises must lead 
to testable predictions, and these predictions must be consistent with experiment and ob- 
servation. Here we will present a theory which appears to satisfy these requirements, but 
which starts with an unfamiliar point of view: There are initially no laws, and instead all 
possibilities are realized with equal probability. The observed laws of Nature are emergent 
phenomena, which result from statistical counting and stochastic fluctuations, together with 
the geography (i.e. specific features) of our particular universe in D dimensions. 



It is likely, of course, that experiment wi 



constraints. For exam 



Higgs-liko bosons |25 



e, supersymmetry 



1 soon confront theory with more stringent 
16 1, SO(d) grand unification jll. [15L 24 ] . and 



27] appear to be unavoidable consequences of the theory presented 



here, but there is as yet no direct evidence for any of these extensions of established physics. 
Many of the ide as p resented here appeared in a more fragmentary and preliminary way 
and conference proceedings [29j. The present paper is self-contained, 



in an earlier paper 



however, and it supersedes all of these precursors. 



II. STATISTICAL ORIGIN OF THE INITIAL ACTION 

Our starting point is a single fundamental system which consists of identical (but dis- 
tinguishable) irreducible objects, which we will call "whits" (because the alternative term 
"bit" has connotations which may be misleading). Each whit can exist in any of N$ states, 
with the number of whits in the ith state represented by rii. An unobservable microstate 
of the fundamental system is specified by the number of whits and the state of each whit. 
An observable macrostate is specified by only the occupancies rii of the states. As discussed 
below, D of the states are used to define D spacetime coordinates x , and Np of the states 
are used to define fields <pk- In this picture, Nature consists of all possible macrostates of 
the fundamental system, which means all possible combinations of and tim, which are 
equivalent to all possible values of the function (tim), or all possible values of <pk (x M ) up 
to phase factors (with the usual convention that x or um here represents the full set of D 
coordinates or occupancies). 
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Let us begin with the coordinates: 

x M = ±n M a (2.1) 

where M = 0,1,. ..,D — 1. It is convenient to include a fundamental length a in this 
definition, so that we can later express the coordinates in conventional units. One can think 
of a as being some multiple of the Planck length t P . 

Notice that positive and negative coordinates correspond to the same occupancies, and 
that this definition is nevertheless physically acceptable: First, two points with coordinates 
that differ by minus signs are typically separated by cosmologically large distances in external 
spacetime. Second, and more importantly, the metric tensor and other physical fields defined 
below need not return to their original values when they are followed along a classical 
trajectory from a point with positive x M to its counterpart with negative x M . I. e., the 
geometry (and other field configurations) described by the two points can be regarded as 
distinct, and in this sense the points themselves are distinct. These different branches of the 
classical geometry (plus fields) are analogous to the branches of a multivalued function like 
z 1 ! 2 or logz, which are taken to correspond to different Riemann sheets, with the points on 
different sheets regarded as distinct. 

Now define a set of fields <p k by 

<Pl(x)= Pk (x) , k = l,2,...,N F (2.2) 

where 

p k (x) = n k (x)/a% (2.3) 

(and x represents all the coordinates). These initial bosonic fields (pk are then real, and 
defined only up to a phase factor ±1. 

Let S (x) be the entropy at a fixed point x, as defined by S (x) = log W (x) (in the units 
used throughout this paper, with ks = h = c = 1). Here W (x) is the total number of 
microstates for fixed occupation numbers n«: W (x) = N (x)!/IL n-i (x)\, with 

N(x) = Y t n i (x) , i = 1,2,..., 7V 5 . (2.4) 

i 

The total number of available microstates for all points x is W = U x W (x) , so the total 
entropy is 

S = J2S(x) , S(x) =logr(jV(aO + l)- J^logTMaO + l) ■ ( 2 - 5 ) 

x i 
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We will see below that n k (x) can be approximately treated as a continuous variable when 
it is extremely large, with 

(99 

— — - = V (N (x) + 1) - V K (z) + 1) (2.6) 
on k (x) 

(> ~ S 1> (1) (N (x) + 1) - ^ (n k (x) + 1) 6 k , k . (2.7) 



<9n fc / (x) dn k (x) 

The functions *0 (z) = dlogT (z) jdz and tp^ 1 ' (z) = d 2 logT(z) /dz 2 have the asymptotic 
expansions 



, ( ^ logz -L.±J^ , ,a, (2H i + _L + £^ ( , 8) 

z=i z=i 
as 2 — )■ oo. It will be assumed that each n k (x) has some characteristic value n k (x) which is 
vastly larger than nearby values: 

n k (x) = n k (x) + An k (x) , n k (x) »> \ An k (x)| (2.9) 

where ">>>" means "is vastly greater than", as in 10 200 >>> 1. This assumption is 
consistent with the fact that the initial action of f!2.27j) and (I2.28P has no lower bound, so 
that n k (x) oo before the extra stochastic term involving (I2.3ip is added. Then it is an 
extremely good approximation to use the asymptotic formulas above and write 

a k (x) An k (x) - ^ a 'k ( x ) l An k (x)] 2 + a '^' ^ Arik ^ Anfc ' ^ ( 2,1 °) 

a k (x) = logiV(x) - \ogn k (x) (2.11) 
a' k (x) = (2n k (x))- 1 -(2N(x)y 1 , a' kk , (x) = (2N (x))" 1 (2.12) 

where N (x) is the value of N (x) when n k (x) = n k (x) for all k, and the higher-order terms 
have been separately neglected in a k (x) and a! k (x). For simplicity, we will also neglect 
the terms involving (2N (x)) . Since there is initially no distinction between the fields 
labeled by k, it is reasonable to assume that they all have the same n k (x) = n(x). It is 
also reasonable to assume that n (x) is nearly independent of x (except possibly in extreme 
situations like the very early universe): n (x) =n and N (x) = N, so that 

o fc (z) = a = log (N/n) (2.13) 
a' k (x) = a' = {2nY 1 . (2.14) 



It is not conventional or convenient to deal with An k (x) and [An k (x)] 2 , so let us instead 
write S in terms of the fields <pk and their derivatives d<j)k/dx M via the following procedure: 
First, we can switch to a new set of points x, defined to be the corners of the Z)-dimensional 
hypercubes centered on the original points x. It is easy to see that 

x . k x . k 

where (■ • • ) in the present context indicates an average over the 2 D boxes labeled by x 
which have the common corner x. Second, at each point x we can write Ank = Ap^a® = 
{{Ap k ) + 5p k ) aft, with (5p k ) = 0: 

" V-A-V-'/.-/ 1 ■ ■ "/'/.- 1 / \ ( ' l>] 

x,k x,k 

S + ^a(Ap fc )a^-^a'[(A Pfc ) 2 + <( ( 5p fe ) 2 >] (a D ) 2 . (2.17) 



S = So + Y, a (( A Pk) + SPk)a°-J2 a '(tt A P^ +6 P^)( a o) 2 ( 2 - 16 ) 



a y\L±p k ) t- \ \op k) ;j \c Ds ~ 

x,k x,k 



Each of the points x surrounding x is displaced by 5x = ±ao/2 along each of the x axes, 
so 

<(« 2 > = (mT) (2.18) 



Ai 



V&; M 2,9(x M ) 2 

E(^" + @) 2 (^) 2 + ^(^) 2 )") (2.20) 

to lowest order, where it is now assumed that at normal energies the fields are slowly varying 
over the extremely small distance Oq. This assumption is justified by the prior assumption 
that n is extremely large: (x) = pk (x) = n& (x) faff implies that 28(fik/(fik ~ Snk/nk and 
4>k = %^ a o D ^ 2 ) so that 5(pk ~ 5rik n k 1 ^ 2 a D ^ 2 . The minimum change in <pk is given by 
onfc = 1 : 50™ in ~ % "o > which means that <50™ m is extremely small if rik is extremely 
large. For example, D = 4, a 1 ~ 10 15 TeV, and n fc ~ 10 200 would give 6<f>f n ~ 1(T 70 
TeV 2 ~ 1(T 38 K 2 . 

In other words, the fields 4>k have effectively continuous values as n — > oo, just as the 
spacetime coordinates have effectively continuous values as ao — > 0. 

For extremely large n it is an extremely good approximation to neglect the middle term 
in (I2.20p . and to replace 2 by 

2 = -p = n/a° (2.21) 



giving 



M 



\dx M J a ° + \d(x M < 2 
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(2.22) 



It is similarly an extremely good approximation to neglect the term in (I2.17P involving 
a 1 (a®) 2 (Ap k ) 2 = (An k ) 2 /2n in comparison to that involving (Apk) clq = (Arik), so that 



where 



2m 2 , 



dx 



M 



d 2 4> k 
16 l<9(x M ) 2 



(2.23) 



m = a 1 , p = m log (iV/n) . (2.24) 
The philosophy behind the above treatment is simple: We essentially wish to replace (f 2 ) 



by (df /dx) 2 , and this can be accomplished because 



f) ~ (ff = ((Sf) 2 ) « <(9//9x) 2 (6x) 2 ) = (df/dx) 2 (a /2f 



(2.25) 



The form of ( I2.23|) also has a simple interpretation: The entropy S increases with the number 
of whits, but decreases when the whits are not uniformly distributed. 
In the continuum limit, 



E 



(12.231) becomes 



s = s + i d u x Y { — 



— Y 

M 



d°X 



(2.26) 



\dx M ) 



16 \d(x M ) 2 



(2.27) 



A physical configuration of all the fields (f>k (x) corresponds to a specification of all the 
density variations Ap^ (x). In the present picture, the probability of such a configuration 
is proportional to W = e s . In a Euclidean path integral, the probability is proportional to 
e~ 5b , where Sh is the Euclidean action for these bosonic fields. We conclude that 



Si, 



-S + constant 



(2.28) 



and we will choose the constant to equal So- 

In the following it will be convenient to write the action in terms of 



,. For 



simplicity, we assume that the number of relevant <fik is even, so that we can group these 
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real fields in pairs to form Nf complex fields $?b,k- It is also convenient to subtract out the 
enormous contribution of <fi by defining 



where ^b is the vector with components *f>b,k and ^b is similarly defined with <pk — > 
the action can be written 



Sb 



d D x 



2mr 



&&1 d^ b a 2 d 2 ^\ 



dx M dx M 



(2.29) 
Then 

(2.30) 



16 d (x M ) 2 d {x M f 

since ^ is constant, with summation now implied over repeated indices like M. 

At this point it is necessary to introduce another fundamental assumption, without which 
the action would have no lower bound: We assume that the whits are randomly perturbed 
by an unspecified environment whose effect can be modeled by an extra term iV ^\^b- Here 
V is a potential which has a Gaussian distribution, with 



V 







V(x)V(x')) = b5(x-x') 



(2.31) 



where b is a constant. Then the complete action has the form 



S, 



d D x 



2mn 



d^ b d% a 2 



a 2 d 2 ^\ d 2 V b 



dx M dx M ' 16 d {x M f d (x M ) 2 



(2.32) 



In the following we will assume that the only fields which make an appreciable contribution 
in (I2.32p are those for which J d D x ^l^b = 0. This assumption is consistent with the fact 
that is constant with respect to all the coordinates and, in the present picture, fields 
corresponding to physical gauge representations have nonzero angular momenta in the 
internal space of Section [V] and Appendices [A] and [B] Then f)2.32p simplifies to 



d D x 



2m 



Q X M Q X M + 



d 2 ^ t 



16 



d (x 1 



-HoVlV b + iV$lV b \ .(2.33) 



III. PRIMITIVE SUPERSYMMETRY 



If F is any functional of the fundamental fields its average value is given by 



' fV$lVV b F 



-Sb [*l* b 



(3.1) 



where (■ • • ) now represents an average over the perturbing potential iV and J T> ^\ V> ^ b is to 
be interpreted as Y[ x k /-^o ^R- e ^ b > k ( x ) I^oo ^ m ^b,k ( x )- The transition from the original 
summation over (x) to this Euclidean path integral has the form (with An = 1 here) 



00 POO POO POO 

V / (n) An / fdn ^ fd (a^ 2 ) 20a? / / 

„_n JO JO Jo 



n=0 

where 0' 



/ #' (3.2) 



m, 



-1/2- 



>, since d ( 



is an extremely good approximation for 



physically relevant fields, and since 0' effectively ranges from —00 to +00. Each then 
becomes a R,e^b,k {%) or Im*^ (x), and the constant factors cancel in the numerator and 
denominator of ( 13. ip . 

The presence of the denominator makes it difficult to perform the average of (I3.ip . but 
there is a trick for removing the bosonic degrees of freedom }fr_ b in the denominator and 



31 



33|: After inte- 



replacing them with fermionic degrees of freedom in the numerator 
gration by parts, f !2.33|) can be written in the form Sb = J d D x \E^v4\I/b. Then, since 



M^]=C(det A) 



-1 



Z>*+D* /e - & [*/'*'] =det A 
where the matrix A corresponds to the operator A and C is a constant, it follows that 

where and ^ r have been combined into 



(3-3) 
(3.4) 

t 

(3.5) 
(3.6) 



and 



S bf = / d D x 



2mn 



dm 




(3.7) 



+ 



<W d 2 ^ 



dx M dx M lQd{x M ) 2 d{ 



x 



(3.8) 

In (13. 7p . tyf consists of Grassmann variables ^f t k, just as consists of ordinary variables 
^b,k, and J V ^ V ^ is to be interpreted as 



dRe^b,k {x 



dim* 



b.k 



d^* fk (x) I dVf tk (x) 



(3.9) 



Recall that and *f?f each have Nf components. 

For a Gaussian random variable v whose mean is zero, the result 



( e -™)=e-^ 2 > 



implies that 



\bjd D x [^{x)^{x)\ 



It follows that 

with 

S E = J d D x 
A special case is 



(F) 



1 

C 



(3.10) 

(3.11) 
(3.12) 

(3.13) 



2m n 



dx M dx M 



+ 



,2 d 2 tft d 2 ^ 



16 d 



1 



d(x M Y 



-/i ^ + -6(^) 2 \ . (3.14) 



(3.15) 



(3.16) 



but according to (HUT) Z = 1, so C = / V^V^/e' SE and 

. . [V&W Fe~ SE 

(F) = 

w fVWVW e~ s ^ 

Notice that the fermionic variables tyf represent true degrees of freedom, and that they 
originate from the bosonic variables ^_ b . The coupling between the fields and ^/ (or 
\j/ b ) is due to the random perturbing potential iV. In the replacement of (13. ip by (I3.16p . F 
essentially serves as a test functional. The meaning of this replacement is that the action 
(I3.14p . with both bosonic and fermionic fields, must be used instead of the original action 
(I2.33p . with only bosonic fields, in treating all physical quantities and processes, if the 
average over random fluctuations in (13.11) is to disappear from the theory. 

Ordinarily we can let ao — > in (I3.14p . so that 



St 



d D x 



(3.17) 



However, the higher- derivative term in the action may be relevant in the internal space 
defined below, where the length scales can be comparable to ao- Also, one should remem- 
ber that the discreteness of spacetime in the present theory implies that a automatically 
provides an ultimate ultraviolet cutoff. 



IV. LORENTZ INVARIANCE 



As mentioned in the first section, within the present theory the laws of Nature arise from 
(1) the statistical counting which leads to ( I2.30p . (2) the stochastic fluctuations of (I2.3ip . 
and (3) the geography (or specific features) of our universe, to which we now turn. 

The most central assumption is that 



(4.1) 



where \l/o is the order parameter for a primordial bosonic condensate which forms in the very 
early universe, and ^>' b represents all the other bosonic fields. The treatment of Appendix |A] 
implies that 







(4.2) 



everywhere. (See f ]Allj) and the comments above (1A2I) and (1A9I) .) The action can then be 
written as 



S E 

Scond 

s b 

Sf 

Sint 

V 



Scond ~\~ Sb + Sf + S. 



>int 



d D x 


1 




2m 


d D x 


1 




2m 


d D x 


1 




2mo 



^—d M ^fd M ^f + (V - /i ) + ^ b (*/*/ 



d D xb { ty f ty f 



(4.3) 
(4.4) 

(4.5) 

(4.6) 

(4.7) 
(4.8) 



In the following, the last term will be neglected in ( 14. 5 p and (|4.6p ; we are thus considering 
the theory prior to formation of further condensates beyond the primordial \EV 

1/2 I 

For a static condensate we could write \l>o = n o Vo, where 7] is constant, rj^rjo = 1, and 
n = \1/q\I/o is the condensate density. This picture is too simplistic, however, since the order 
parameter can exhibit rotations that are analogous to the rotations in the complex plane of 
the order parameter ip s = e %es nj 2 for an ordinary superfluid: 



#o (x) = U (x) n (x) 1/2 



r]o 



ulu. 



1 . 



(4.9) 
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After an integration by parts in (14. 4p (with boundary terms always neglected in the present 
paper), extremalization of the action gives the classical equation of motion for the order 
parameter: 

1 -d M d M % + (V - fM ) % = . (4.10) 



2m 

(In this equation the ordinary field densities *ff' b ty' b and ^\^f are neglected in comparison 
to the primordial condensate density ^q^o-) 

In specifying the geography of our universe, it will be further assumed that \l/o can be 
written as the product of a 2-component external order parameter ty ex t, which is a function 
of 4 external coordinates x^, and an internal order parameter \l/; n t, which is primarily a 
function of d internal coordinates x m , but which also varies with x M : 

$ = ^(X") V mt (x m ,xn (4.11) 

* M = U ext (x») n ext (^) 1/2 Vext , n = 0, 1, 2, 3 (4.12) 
^ int (x m ,x^) = U int (x m ,x^n int (x m ,x^) 1/2 Vint , m = 4,...,D-l (4.13) 

where again r] ext and r] int are constant, and T)l x tVext = VintVint = 1- Here, according to a 
standard notation, x M actually represents the set of x M , and x m the set of x m . 
Let us define external and internal "superfluid velocities" by 

mov^ = -iU~]d^U ext , m v m = -iU^]d m U int . (4.14) 

The fact that Ul xt U ext = 1 implies that {d^Ul x ^ U ext = —U\ xt (p^Uext) with U\ xt = U~ xt , or 
m v^ = i (dpUlgt) U ext , so that 

vl = v,. (4.15) 
For simplicity, let us first consider the case 

d»U int = (4.16) 

for which there are separate external and internal equations of motion: 

with 

Vint = 1*0 ~ ^ext ■ (4.18) 

The quantities fi int and V have a relatively slow parametric dependence on x^. 
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When fl4~T2l) and (gjg) are used in (Q7j) . we obtain 



t V 2 



nli 2 ^ = (4.19) 



and its Hermit ian conjugate 



t 1/2 



2 moV " - 2^ 



l/ 2 n 



Subtraction gives the equation of continuity 



9^ = 



and addition gives the Bernoulli equation 

1 ,2 



3a ~ n ext Vext v pJ]ext 



2 m V ext + P ext = Vext 



where 



ext 



Vext V n V n Vext 



ext 



2m ( 



1/2 



(4.20) 



(4.21) 



(4.22) 



(4.23) 



Since the order parameter ^ e xt in external spacetime has 2 components, its "superffuid 
velocity" can be written in terms of the identity matrix cr° and Pauli matrices o a : 



(4.24) 



//, a = 0, 1, 2, 3 . 

It will be assumed that the coordinate system can be chosen such that 

v° k = v% = , fc,a = l,2,3 

to a good approximation, yielding the simplification 

-m v%v* + P ext = fx ext . 

As v" varies, \i ext varies in response, with fi int determined by ( 14. 18ft . 

Now expand ty' b in terms of a complete set of basis functions ipl nt in the internal space: 



(4.25) 



(4.26) 



V b W W nt (X" 1 ) 



with 



2mn 



d m d m - ii int + v (* m ) = ^ (* m ) 



(4.27) 

(4.28) 
(4.29) 
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and with the usual summation over repeated indices. For reasons that will become fully 
apparent below, but which are already suggested by the form of the order parameter, each 
4>l {x 11 ) has two components. As usual, only the zero (e r = 0) modes will be kept. (To 
simplify the presentation, the higher-derivative terms are not explicitly shown in the present 
section; they will be restored in the next section.) When f l4.27p - fl4.29p are then used in (14. 5 p 
(with the last term neglected), the result is 



S b = I d 4 x 4 ( -^fyfy - Vext ) ^ (4.30) 



where ip b is the vector with components ip r b . 
Let ipb be written in the form 

A(x») = U ext (x»)^ b (xn (4.31) 

or equivalently 

V b {x») = U ext {x»)V b W • (4.32) 

Here ipt, has a simple interpretation: It is the field seen by an observer in the frame of 
reference that is moving with the condensate. In the present theory, a (very high density) 
condensate \I/o forms in the very early universe, and the other bosonic and fermionic fields are 
subsequently born into it. It is therefore natural to define the fields tpl in the condensate's 
frame of reference. 

Equation (I4.3ip is, in fact, exactly analogous to rewriting the wavefunction of a particle 
in an ordinary superfluid moving with velocity v s : ip par (x) = exp {imv s x) ip par (x) . Here ip par 
is the wavefunction in the superfluid's frame of reference. 

When (147311 is substituted into (147301) . the result is 

If n s and are slowly varying, so that P ext and <9 M t> M can be neglected, (14.251) and ( 14.261) 
lead to the simplification 

S b = - j d 4 x 4 C^9& + iv^dA ij b . (4.34) 

In the remainder of the paper it will also be assumed that 

P„ < rnov* (4.35) 
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S b = d A x ipl 



^ b ■ (4.33) 



for states with momenta p u corresponding to energies ~ 1 TeV or less (as would be the case 
if we had, e.g., m = % 1 ~ 10 15 TeV and v£ ~ 1), so that ( I4.34p reduces to just 

S b = J d 4 x 4ie^a a d^ b , e£ = . (4.36) 

Up to this point there has been no distinction between contravariant and covariant vectors, 
and could have been written, e.g., t>^, but it is important beyond this point that be 
treated clS 8u contravariant 4- vector. 

The above arguments also hold for fermions, with 

S f = J d D x ^--L-tf f d M d M * f - fio*\*f + W}tf/J (4.37) 

V f (x»,x m )=r f (xnV int (x m ) (4.38) 

leading to the final result 

S f = J d A x ^ie^d^f . (4.39) 

The present theory thus yields the standard action for Weyl fermions, with the gravi- 
tational vierbein interpreted as essentially the "superfluid velocity" associated with the 
condensate \I/o- The path integral still has a Euclidean form, and the action for bosons is 
also not yet in standard form, but we will return to these points below. 



V. GAUGE FIELDS 



Let us now relax assumption (14.161) and allow Ui nt to vary with the external coordinates 
x^. Equation (14.101) is satisfied if (14.171) is generalized to 



—d^ - fi ext ) V ext (a;") ^ 
2m / 



int (•£ ) 3^^) 



with ^i n t required to satisfy the internal equation of motion (at each x A 



(5.1) 



E— 



d 2 



+ 



a 2 d 4 



V (x m )-fi mt 







(5.2) 



8(x m ) 2 ' I6d(x m ) 4 

The higher-derivative term of (I3.14p has been retained and two integrations by parts have 
been performed. (In order to simplify the notation, we do not explicitly show the weak 
parametric dependence of /ij nt , Vo, and rii nt on x M .) This is a nonlinear equation because (at 
each x^) V (x m ) is mainly determined by n int = ^>\ nt ^> int . 
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The internal basis functions satisfy the more general version of (I4.28P with e r = 0: 



^(* m ,^) = 0. (5.3) 



This is a linear equation because Vo (a: m ) is now regarded as a known function. 

The full path integral involving ( 13.14)) contains all configurations of the fields, including 
those with nontrivial topologies. In the present theory, the geography of our universe includes 
a topological defect in the d- dimensional internal space which is analogous to a vortex. (See 
Appendix |Al) The standard features of four- dimensional physics arise from the presence of 
this internal topological defect. For example, it compels the initial gauge symmetry to be 
SO{d). 

The behavior of the condensate and basis functions in the internal space is discussed in 
Appendices lAl and iBl In ( 1A15j) . the parameters (p i specify a rotation of \I/j ni (x m ,x M ) as the 
external coordinates x M are varied, and according to (IA16[) the Jj satisfy the SO(d) algebra 

Jilj - JjJi = iCip k . (5.4) 

For simplicity of notation, let 

(r \Q\r') = J d d x$lQ$ nt with (r \r') = 5 rr> (5.5) 
for any operator Q, and in particular let 

tf = {r\li\r') (5.6) 
with the matrices t r {' (which are constant according to (1A17I) ) inheriting the SO(d) algebra: 

fa - t 3 t t y r ' = j2( r 1 7 y) ( r "\ w) - E ( r i 7 >"> < r "i J y) ^ 



= {r\JiJj\r')-{r\J j J i \r') (5.8) 

= idtfk ■ ( 5 - 9 ) 

The ti are the generators in the A^-dimensional reducible representation determined by 
the physically significant solutions to (I5.3p . which spans all the irreducible (physical) gauge 
representations. 

When x M — > x M + dx^, ^ int and ip r int rotate together, and f)A15j) implies that 

d^\ nt [x m ^) = 4^^^,/) (5.10) 



= -i^J^L^,^) (5.H) 
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where 



A] 



(5.12) 



The A 1 ^ will be interpreted below as gauge potentials. In other words, the gauge potentials 
are simply the rates at which the internal order parameter ty int (x m , x^) is rotating as a 
function of the external coordinates x^. 

Let us return to the fermionic action ( 14. 37ft . If ( 14. 38ft is written in the more general form 



9, (x», x m ) = r f m (* m , = u ext (xn r f m r mt (* m , ^) 



we have 



d^ f = U ext (a?) (9; - im e»a a - i A\J, L ) Vi 
where the prime indicates that does not operate on ipl nt , and 

d d x¥ f d^ f 



r 

int 



(5.13) 



(5.14) 



d d x ^1 K - im e^ a - i A^J t ) - im e»,a<*' - i Aj J, ) « t (5.15) 
^ (r | (9; - mo e^ a - i Apt ) ]T |r") (r"| - mo e> a ' - * < 7* ) |r') < (5.16) 



5 rr // (<9 M - im e£cr a ) - iA 1 ^" \ |<W ( ^ - im e£,er 



Then (I4.37[) becomes 



< (5.17) 
(5.18) 



where 



D^d^-iA^. 



With ( 14.251) and the approximations above (14.341) . ( 14.261) implies that 



5/ 



d^xip) 



(5.19) 



(5.20) 



7 V 2 W 

This is the generalization of ( I4.34[) or (14.391) when the internal order parameter is permitted 
to vary as a function of the external coordinates x M . Again, for momenta and gauge potentials 
that are small compared to moe^ (as would be the case at normal energies and fields if m,Q 



16 



were comparable to a Planck energy and were ~ 1), the first term may be neglected. 
Furthermore, the entire treatment above can be repeated for the bosonic action, finally 
giving 

S f = J d 4 x^ f ie»a a D^ f , S b = J d^x^ie^D^ . (5.21) 

VI. TRANSFORMATION TO LORENTZIAN PATH INTEGRAL: FERMIONS 

All of the foregoing is within a Euclidean picture, but we will now show that, in the 
case of fermions, there is a relatively trivial transformation to the more familiar Lorentzian 
description. A key point is that the low-energy operator ie^a a D^ in Sf is automatically 
in the correct Lorentzian form, even though the initial path integral is in Euclidean form. 
It is this fact which permits the following transformation to a Lorentzian path integral. 
Within the present theory, neither the fields nor the operators (nor the meaning of the time 
coordinate) need to be modified in performing this transformation. 

In a locally inertial coordinate system, the Hermitian operator within Sf can be diago- 
nalized to give 

S f = J d A x ^}(x) ia^D^jix) (6.1) 
= J2 S ^f(s)a(s)^ f (s) (6.2) 

where 

ijf (x) = J2 u ( x > s ) ^7 ( s ) > = / d ' x c/t ( x > s ) ^ (*) ( 6 - 3 ) 

s J 

with 

i^DpU (x,s) = a(s)U(x,s) (6.4) 
/ d?xW (x,s)U (x,s') = 5 SS , , ^2u(x,s)U*(x',s) = 6(x-x') . (6.5) 

J s 

Here, and in the following, x represents a point in external spacetime, and U (x, s) is a 
multicomponent eigenfunction. There is an implicit inner product in 

C/t ( x> s ) ^ (X) = U r & S) i>f, r (X) (6.6) 

r 

with the 2N g components of ipf (x) labeled by r = 1, N g (spanning all components of all 
irreducible gauge representations) and a = 1,2 (labeling the components of Weyl spinors), 
and with s and (x, r, a) each having N values. 
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Evaluation of the present Euclidean path integral (a Gaussian integral with Grassmann 
variables) is then trivial for fermions; as usual, 

(6.7) 
(6.8) 

(6.9) 

s 

with 



Vipj (x) T>ip f (x) e~ s f 

II / d ^f,ra( X ) / dll)f >ra (x) e~ S f 
x,ra 

II*/00 



z f I s ) 



a{s) 



since the transformation is unitary 



34|. Now let 



where 



z f [s) 



ia (s) 



so that 



c f Z f 



n 



(6.10) 
(6.11) 



(6.12) 
(6.13) 



(6.14) 
(6.15) 

(6.16) 



This result holds for the path integral over an arbitrary time interval, with the fields, oper- 
ator, and meaning of time left unchanged. 

The transition amplitude from an initial state to a final state is equal to the path integral 
between these states, so transition probabilities are the same in the Lorentzian and Euclidean 
descriptions. This result is consistent with the fact that the classical equations of motion are 
also the same, since they follow from extremalization of the same action. Furthermore, using 
the method on pp. 290-291 or 302-303 of Ref. 34j, it is easy to show that the magnitude 
\G (x,x')\ of the 2-point function is again the same, so particles propagate the same way in 
both descriptions. This result is also obtained in Appendix O with a different method. 
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When the inverse transformation from ipj to ipf is performed, we obtain 

Zf = ( Vip) (x) V^P f (x) e iS f (6.17) 



with Sf having its form ( 16. ID in the coordinate representation. 

One may perform calculations in either the path-integral formulation or the equivalent 
canonical formulation, which can now be obtained in the standard way: Let us use the 
notation to indicate that the fields in a path integral are specified to begin in a state \a) 
at time t a and end in state \b) at time t b , and also to indicate that a path integral showing 
these limits has its conventional definition (so that it may differ by a normalization constant 
from Zf as defined above). Then the Hamiltonian Hf is defined by 

(b\U f (t b ,t a )\a) = / Vi,){x) Vif> f (x) e lS ? (6.18) 

J a 

ijU f (t, t a ) = H f (t) U f (t, t a ) , U f (t a , t a ) = 1 (6.19) 



as in (9.14) of Ref. 34j. I.e., the time evolution operator Uf (tb, t a ) is defined to have the same 



effect as the path integral over intermediate states, and it is then straightforward to reverse 



the usual logic which leads from canonical quantization to path-integral quantization 34LI35]. 



VII. TRANSFORMATION TO STANDARD FIELDS AND LORENTZIAN PATH 
INTEGRAL: BOSONS 

For bosons we can again perform the transformation (16.31) to obtain 

s b = Yl ^» a 00 ^ 00 • ( 7 - x ) 

• • s 

The Euclidean path integral is 

Z b = JviPi{x) Vij> b (x) e~ Sb (7.2) 
d (Re^ ra (x)) / d(lm^ ra (x)) e- Sb (7.3) 
= l[z b (s) (7.4) 

s 

with 

p oo r oo 

z b (s)= d(Re^ b (s)) d(lm^ b (s))e~ Sb . (7.5) 



oo 
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We will now show how this action can be put into a form which corresponds to scalar 
bosonic fields plus their auxiliary fields. First, if the gauge potentials A 1 ^ were zero, we would 
have 

ta^d^U (x, s) = a (s) U° (x, s) . (7.6) 

Then 

U°(x,s)=V' 1/2 u(s)e ip ^ , p s -x = n^rix" , 7^ = diag(-l,l,l,l) (7.7) 

(with V a four-dimensional normalization volume) gives 

- v^P'sU (x, s) = a (s) U° (x, s) (7.8) 

where <r M implicitly multiplies the identity matrix for the multicomponent function U° (x, s). 
A given 2-component spinor u r (s) has two eigenstates of p k s a k \ 

p k s a k ut (s) = \f s \ut(s) , p k a k u; (s) = - \f s \ u~ (s) (7.9) 

where p s is the 3- momentum, with magnitude The multicomponent eigenstates of 

ia^d^ and their eigenvalues a (s) = p°=F \~P s\ thus come in pairs, corresponding to opposite 
helicities. 

For nonzero A 1 , the eigenvalues a (s) will also come in pairs, with one growing out of 
do (s) and the other out of its partner a (s r ) as the A 1 ^ are turned on. To see this, first write 
flEU as 

(id + A&i) U (x, s) + a k (id k + A\ti) U (x, s) = a (s) U (x, s) (7. 10) 



or 



id 5 rr i + A\tY' j U r i (x, s) + P rr iU r i (x, s) — a (s) 8 rr 'U r i (x, s) = (7-11) 
P rr , = a k (id k 8 rr , + A{tf) (7.12) 



with the usual implied summations over repeated indices. At fixed r, r' (and x,s), apply a 
matrix s which will diagonalize the 2x2 matrix P rr i, bringing it into the form p rr >a 3 +p rr ,a , 
where p rr i and p rr , are 1-component operators, while at the same time rotating the 2- 
component spinor U r r. 

sP rr ,s _1 = P'„, = p rr >a 3 +p rr ,a° , U' r , = sU r > (7.13) 

--'::)• 
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But P rr > is traceless, and the trace is invariant under a similarity transformation, so p rr , = 0. 
Then the second term in (17. lip (for fixed r and r') becomes s~ 1 p rr >a 3 U^, (x,s). The two 
independent choices 

U',(x,s) oc [ 1 I , a 3 U' r , (x, s) = UL (x, s) (7.15) 

W 

U' r , (x,s) oc , a 3 U', (x, s) = -U' r , (x, s) (7.16) 

give ±s~ lr p rT iU' r , (x,s). Now use s _1 L^, = U r i to obtain for (17. lip 

(^id 5 rr ' + A\t™' \ U r i (x, s) ± p rr /U r > (x, s) — a (s) 5 rr rU r > (x, s) — (7-17) 

so 117. 10p reduces to two sets of equations with different eigenvalues a (s) and a (s')\ 

a (s) = a,\ (s) + a 2 (s) , a (s') = a\ (s) — a 2 (s) (7-18) 

where these equations define a± (s) and a 2 (s). Notice that letting a k — > —a k in (I7.10p 
reverses the signs in (I7.17p . and results in a (s) — > a (s')\ 

(id + AlU) U (x, s) - a k (id k + U (x, s) = a (s') U (x, s) . (7. 19) 

The action for a single eigenvalue a (s) and its partner a (s') is 

s b (s) = ft (s) a (s) (s) + ft (s') a (s') ft (s') (7.20) 
= ft (s) (a x (s) + a 2 (s)) ^ 6 (s) + ^ (s') (a x (s) - a 2 (s)) ^ 6 (s') . (7.21) 

For ai (s) > 0, let us choose a 2 (s) > and define 

fo(s') = a( S ) 1/2 b ( S ) = (a 1 ( S ) + a 2 ( S )) 1/2 b ( S ) (7.22) 
^ b (s) = a ( S )" 1/2 F b (s) = ( ai (s) + a 2 (s))- 1/2 F 6 (s) (7.23) 



so that 



s 6 (s) = 0;(s)a( S )0 b ( S ) + F;( S )F b (s) , oi(s)>0 (7.24) 



where 



a(s) = a (s) a (s') = a x (s) 2 — a 2 (s) 2 . (7-25) 
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For a\ (s) < 0, let us choose a 2 (s) < and write 

^b(s') = (-a(s)) 1/2 fe (s) = (-oi(s) -a 2 (s)) 1/2 fe (s) 

= (-a( S ))- 1/2 F b ( S ) = (- ai ( S )-a 2 ( S ))- 1/2 F b ( S ) 



so that 

3& (s) = - 
Then we have 

s b = s b (s) 

* ' s 

I 

= [T b (s)a(s)^ b (s)+F* b (s)F b (s) 

oi(s)>0 



6 *(5)a( S )0 6 (5)+F;(5)F 6 ( S ) 



, a a (s) < . 



(7.26) 
(7.27) 

(7.28) 
(7.29) 



[fb (*) a 00 06 00 + 00 ^ 00 (7-30) 

ai(a)<0 



where a prime on a summation or product over s means that only one member of an s, s' 
pair (as defined in (J77T7J) and ( 17.1 8p ) is included. 

Each of the transformations above from ip b to <p b and F b has the form 



A(s') = A(sf 2 <j )b (s) , tp b (s) — A (s)~ l/J F b (s) 



1/2" 



(7.31) 



so that dip b (s') = A (s) d<p b (s) , dip b {s) = A(s) 1 ^ 2 dF b (s), and the Jacobian is 
Y\' s A (s) 1 ^ 2 A (s) -1 ^ 2 = 1. These transformations lead to the formal result 



z b = n^( s ) • n^oo ■ n z f( s )- n ^oo 

oi(s)>0 ai(s)<0 



(7.32) 



s>0 



s<0 



where 



s<0 < — ► a (s) = di (sf - a 2 (s) 2 < if a\ is) > 
^ — )■ a (s) = ai (s) 2 — a2 (s) 2 > if ai (s) < 



(7.33) 
(7.34) 



with s > otherwise, and where 

^( s ) = /°°d(Re0 6 ( S )) /' OO C /(Im0 f) ( S ))e T|S(s)| [( Rc ^ (s) ) 2+ ( Im?bW ) 2 ] (7.35) 
J —00 J —00 

Zf ( s )= f°° d(ReF b (s)) /"°°rf(ImF b ( S ))e T [( Rc ^ (s) ) 2+ ( Im ^ {s) ) 2 ] (7.36) 

J ~oo J —oo 

with the upper sign holding for s > in (I7.35P and a\ (s) > in ( 17.361) . (Also, a% (s) = is 
treated as a limiting case of a± (s) > 0, and a(s) =0 requires the usual convergence factor 
in the path integral.) 
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The above formal expressions would be highly divergent for the modes with the lower 
signs. If represents the field consisting of these modes, it therefore cannot be treated 
with the approximations used above. Its behavior is instead controlled by second-derivative 
and nonlinear terms with the same form as in (14.51) . which respectively restrict the number 
of modes which have negative action and the contribution of each such mode. One might 
conjecture that (^f' b ^f b \ behaves effectively as a charge-neutral condensate density, very 



roughly analogous to the QCD vacuum condensate; but the only certainty is that treating 
the components of this field properly requires methods beyond those used in the present 
paper, so these modes will be omitted in the following discussion. 
Then (I7.32p is reduced to 



Z b = Y[z<f>(s) ■ JJ z F (s) , z$ (s) = , z F (s) = n (7.37) 

s>o oi(s)>o ' a ^ ' 



where s > means 

ai (s) > \a 2 0)| ifai(s)>0 , \ai (s)| < \a 2 (s)| ifoi(s)<0 (7.38) 

according to the definition below ( 17.321) . Recall that if the gauge potentials A 1 ^ were zero, 
we would have a\ — u and a 2 = \~$\, where u is the frequency and p the 3- momentum. 



Now let 



where 



z L As) 



Zl = jv4 (s)VM*) VF\ (s)VF b (s) e iS » (7.39) 

= n n (7.40) 

s>0 ai(s)>0 



d(Re^ b (s)) / OO rf(Im0 6 ( S ))e jS(s) [( Re0&(s) ) 2+ ^ Im0b(s) ) 2 ] (7.41) 



— oo 
7T 



i^j-r (7.42) 
a{s) 

/OO p OO r , \ 2 / \ 2~| 

d(ReF b (s)) / d(ImF b (s))e^ RcFb{s) > + ( lmFb(s) > \ (7.43) 
oo J —oo 

= in (7.44) 

since f^dx f^dy exp ( ia (x 2 + y 2 )) = in /a. (Nuances of Lorentzian path integrals are 
discussed in, e.g., Ref. [34], P- 286.) We have then obtained 

Zt = c b Z b (7.45) 
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where c b is a product of factors of i. 

To return to the coordinate representation, let 



®b (x) = U ( x > s ) fa ( s ) 



<t>b ( s ) 



s>0 



Fh (x) 



U(x,s)F b (s) , F b (s) 

ai(s)>0 



d 4 x U ] (x, s) $ 6 (x) 
d A x W (x, s) T b (x) 



(7.46) 
(7.47) 



Recall that, according to (17.101) and (I7.19p . 

ia^D^U (x, s) = a (s) U (x, s) , ia^D^U (x, s) = a (s') U (x, s) 



(7.48) 



with a = cr°, o k = —<J k , a (s) = a x (s) + a 2 (s), and a (s') = a\ (s) — a 2 (s), so the action in 
(I7.39P becomes (with a (s) = a x (s) 2 — a 2 (s) 2 ) 



s>0 



oi(s)>0 



J d 4 x U>\ (x) io^D^ io»'D^ b (x) + U\ (x) io^D^ ia^D^ b (x) + J* (x) T b (x) 



(7.49) 



.(7.50) 



Now, just as an ordinary vector can be made to align with a single coordinate axis through 
a rotation, Q b (x) can be rotated (at each point x) to have only one component <&' bR (x) in 
each irreducible gauge representation R, with A\ — > A Rk to keep the action invariant. (This 
is similar to a transformation to the unitary gauge of SU(2).) Although <&' bR (x) itself still 
has 2 components, the nonabelian term in (17.501) can now be diagonalized: 

KW^Ltr^ k 'A'^t^' b ] R = h' b % (a k a k 'A'^A^ + a k V A^t?) $^(7.51) 

(7.52) 
(7.53) 



2 



1 Rk°i' ^b,R ■ 



If we neglect d^b {x), it is not difficult to show that the other terms can be similarly treated 
and the two components of are decoupled. After an inverse rotation within (17. 53ft in each 
irreducible representation, ( 17.50P then becomes 



S b 



d*x 



$J (x) n^D^i (x) + P h (x) T h (x) 



(7.54) 



The two components of (x) or T b (x) are not independent, since they are both determined 
by 06 ( s ) or F b (s). 
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The return to the coordinate representation is much simpler when ia^D^ia^' D^i and 
rfD^Pv have the same spectrum of eigenvalues, so that there are one-component orthonor- 
mal eigenfunctions U b (x, s) satisfying 



r]^D^D v U b {x,s) = a{s)U b {x,s) . 



(7.55) 



In this case, let 

4> b (x) 



Y^U b (x,s) (j) b (s) , 4> b (s) = / d*x U ] b (x, s) 4> b (x) 

s>0 



(7.56) 



F b (x) = U » ^ s ) Fb ' Fb = / d ' x U b ( x > s ) F b ( x ) ( 7 - 57 ) 

ai(s)>0 

(with the number of values of x chosen to match the number of values of s for a unitary 
transformation), so that 



Sh = / d 4 x 



ij (x) if v DyP v <$> h (x) + (x) F b (x) 



with 



V d>t (x) V <p b (x) V F ] h (x) V F h far) e iSb . 



(7.58) 



(7.59) 



Again, this is the path integral for an arbitrary time interval, and one can define a time 
evolution operator and Hamiltonian as in Section I VI I Notice, however, that the variations 
in <pb{x) and F b (x) are constrained by the expansions of (17.561) and ( 17. 57(1 . These funda- 
mental scalar bosonic fields should therefore exhibit nonstandard behavior. For example, 
as mentioned at the end of Appendix [O, predictions based on radiative corrections will be 
altered. 

VIII. SUPERSYMMETRY AND GRAVITY 



The total action for fermions and bosons is 



Sf + Sb = I d A x 



^ f (x) ia^D^f (x) + 4 (x) v^D,D u <p b (x) + f£ (x) F b (x)] (8.1) 



which in a general coordinate system becomes 



Sf + Sb 



d x e 



(x) ie£ a a D^ i> (x) - sT [D^ (x) ) D u 4> (x) + F^ (x) F (x) 



(8.2) 
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where is the metric tensor and 

e = det el = (- det g^) 1 ' 2 , , D„ = D^ + e~ l/2 d^e 1/2 (8.3) 
^ ( x ) = e ~ 1/2 ifj f (x) , (j)(x)= e" 1/2 <j) b (x) , F (x) = e~ 1/2 F b (x) . (8.4) 

We thus obtain the basic form for a supersymmetric action, where the fields 0, F, and 
ip respectively consist of 1-component complex scalar bosonic fields, 1-component complex 
scalar auxiliary fields, and 2-component spin 1/2 fermionic fields. These fields span the 
various physical representations of the fundamental gauge group, which must be SO(d) in 
the present theory, where d = D — 4. (More precisely, the group is Spin(d), but SO(d) is 
conventional terminology.) I.e., tp includes the Standard Model fermions and Higgsinos, and 
<p corresponds to the sfermions and Higgses. 

According to ( 18. 2ft . the coupling of matter to gravity is very nearly the same as in standard 
general relativity. However, if the action is written in terms of the original fields ipf and 4>b, 
there is no factor of e. In other words, in the present theory the original action has the form 
J d A x C, whereas in standard physics it has the form J d A x e C. For an C corresponding to 
a fixed vacuum energy density, there is then no coupling to gravity in the present theory, 
and the usual cosmological constant vanishes. This point was already made in an earlier 



paper 28(], where the "cosmological constant" was defined to be the usual contribution to 
the stress-energy tensor from a constant vacuum Lagrangian density Co, which results from 
the factor of e. However, as was also pointed out in this 1996 paper, "There may be a much 
weaker term involving SCo/dg^, but this appears to be consistent with observation." 
In the present theory, the Einstein-Hilbert action, with 



C G = \r p 2 e^R, (8.5) 



and Maxwell- Yang-Mills action, with 



C 9 = -\g 2 eF; u F; a g^g™ , (8.6) 



must arise from the response of the vacuum. (Here go is the coupling constant for the 
fundamental gauge group and £ 2 P = 8ttG.) I.e., these terms are interpreted as essentially 



arising through a Sakharov "induced gravity" mechanism 36l-l41|. The gravitational and 
gauge curvatures must ultimately originate from 4-dimensional topological defects associated 
with the vierbein of (I4.36P and gauge potentials of (I5.12p . 
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IX. CONCLUSION 



Although a complete theory will require much additional work - including more detailed 
predictions for experiment - the following have been shown to arise as emergent properties 
from the initial statistical picture: Lorentz invariance, the general form of Standard-Model 
physics, an SO(d) fundamental gauge theory (with e.g. SO (10) permitting coupling con- 
stant unification and neutrino masses), supersymmetry, a gravitational metric with the form 
(—,+,+, +), the correct coupling of matter fields to gravity, vanishing of the usual cosmo- 
logical constant, and a mechanism for the origin of spacetime and quantum fields. 

The new predictions of the present theory appear to be subtle, but include Lorentz 
violation at very high energies and nonstandard behavior of fundamental scalar bosons. 

Appendix A: The internal space 

The internal space of Section |V] is d- dimensional, with an SO(d) (or more precisely 
Spin(d)) rotation group and its vector, spinor, etc. representations - for example, the 10 
and 16 representations when d = 10. It may be helpful to begin with an analogy, however, 
in which external spacetime is replaced by the z-axis. The internal space is replaced by 
an xj/-plane, with internal states described by 2-dimensional vector fields (rather than the 
higher-dimension vector and spinor fields considered below). One of these states is occupied 
by the condensate, and is represented by a vector V\ which points radially outward from 
the origin at all points in the xy-plane when z — 0. The other state is an additional basis 
function, represented by a vector v 2 which is everywhere perpendicular to V\. But v% is 
allowed to rotate as a function of z, so it has both radial and tangential components after a 
displacement along the z-axis. Then v 2 is forced to rotate with V\ - i.e., the condensate - 
in order to preserve orthogonality. 

Now let us turn to the actual internal space, first considering a set of d- dimensional vector 
fields ip r vec . Let ip® ec represent the state occupied by a bosonic condensate. In the simplest 
picture, and at some fixed Xq, only the rth component of the field ip r vec is nonzero along 
some radial direction in the internal space, making the fields trivially orthogonal in that 
direction. Then, with still fixed, ip r vec (x m ) in all other radial directions is obtained from 
the original ip r vec (x™) by rotating it to x m . In other words, the field at each point in the 
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internal space is identical to the field that would be obtained at that point if the original 
field ipl ec (x™) were subjected to a rotation about the origin. This produces an isotropic 
configuration for the condensate and each basis function. As in (I4.13P we can write 

r vec (x m ) = U vec (x m , O ^ ec «) . (Al) 

Just as in the analogy, a field that is radial at x™ will also be radial at all other points x m . 
However, a general tpl ec (x™) permits a general vortex-like configuration of the condensate. 

Also as in the analogy, the state ip^ec °f the condensate is allowed to rotate as a function 
of x M (because such a rotation does not alter the internal action). Since the other basis 
functions ipvec are required to remain orthogonal to ip® ec and each other, they are required 
to rotate with the condensate. Then (IA1I) becomes more generally 

^vec ( X ' xfl ) = U vec (x , X ; X^, Xq) Ipvec ( X ) X o) (^2) 

with 

rj c (x m , x*) v vec (x m , a?) = rj c w, 4) v vec «, o = s rr , (as) 

since 

^vec i X i X i X ^ i ^o) Uvec { x j ^0 i X ^ i X o) = ^ ■ C^^) 

In general (with x M fixed), let ip (x) represent a multicomponent basis function with 
angular momentum j at a point x in the d-dimensional internal space. After a rotation 
about the origin specified by the d x d matrix R, it is transformed to 

i)' (x) =K(R)$ (R^x) (A5) 

where 1Z (R) belongs to the appropriate representation of the group Spin (d). However, we 
require that the field be isotropic, so that it is left unchanged after a rotation: 

$(x) =i>{x) . (A6) 

Then we can define ip (x) at each value of the radial coordinate r by starting with a ip i x o) 
and requiring that 

ip{x) =K(R) f(x ) , x = Rx . (A7) 
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With this definition, tp (x) is a single- valued function of the coordinates only if j is an integer. 
If j = 1/2, e.g., ip (x) acquires a minus sign after a rotation of 2n, but it is single-valued on 
the Spin (d) group manifold. 

Multivalued functions are well-known in other similar contexts, such as the behavior of 
the phase of an ordinary superfluid order parameter ip s = e %es nj 2 around a vortex, which 
becomes discontinuous if it is required to be a single- valued function of the coordinates 4l| . 
In the same way, z 1 ^ 2 exhibits a discontinuity across a branch cut if it is required to be a 
single- valued function and z is restricted to a single complex plane. I.e., 

z l/2 = \ z \V* e i*/2 

1 /2 1 /2 

gives + \z\ ' for = and — \z\ ' for <\> = 2ty. But when defined on a pair of Riemann 
sheets, z x l 2 is a continuous function, and the same is true of ip (x) as we have defined it 
above, on the group manifold. The key idea in either case is to extend the manifold over 
which the function is defined, so that there are no artificial discontinuities. 

A vectorial condensate and vectorial basis functions are appropriate for the simplest 
Higgs-like fields and their superpartners. Similarly, spinorial fields tp r are appropriate for 
ordinary fermions, sfermions, and a possible primordial condensate occupying a state ip%- 
(In the present context, of course, "vector" and "spinor" refer only to properties in the 
internal space.) Again, let ip r sp (a;™) represent a field along some radial direction in the 
internal space at some fixed Xq. Then the field configuration for every point x m is obtained 
by taking ip r sp (x m ) to be identical to the field that would be obtained at that point if ip r sp {x™) 
were subjected to a rotation, with 

ti p (x m ) = U sp (x m ,x™)^ p (x™) (A8) 

as in fTA~7j) . 

Again, the state ifjg p of the condensate is allowed to rotate as a function of x M , and since 
the other basis functions ipl p must remain orthogonal to %fj° sp they are required to rotate with 
the condensate. The general version of ( 1A8I) is then 

^Sp ( X l Xfl ) = U sp (x , Xq ; Xq) ljj gp {Xq , Xq) . (^9) 

The same reasoning applies to each irreducible representation, and thus to the combined 
set of fields i\f int (x m ,x M ): 

(*'"\ = U mt (x' m , x m ; x">, x») C t x») (A10) 
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with 

€1 tint *"*) = tilt tint (*"\ *") = Srr> • (All) 

So that the internal action will be unaffected as x^ — > x ,fl , we require that the order 
parameter experience a uniform rotation, described by a matrix TZint which is independent 
of x m . Then Umt has the form 

U mt (x ,m , x m ; x"\ x^) = U mt (x'», a?) U mt (x ,m , x m ) . (A12) 

(Notice that (1A12j) is to be distinguished from a rotation about the origin, which is given 
by flA5p . and which according to ( ]A6j) would leave ip (x m ) unchanged rather than rotated at 
each point x m .) It follows that 

'Vint ( x ™ \ = Hint (x^, Xq) ^)\ nt (x™ , X ) . (A13) 

We define the parameters 5</>j by 

K int (x" + 5x^ xft) = TZ int (x», xft) (l - z <ty, Ji) (A14) 

or 

^L^" 1 ) = -iHi'Ji'ilntix™) as a^-^ + fo" (A15) 

Jj = lZi n t (x^ , Xq) JilZ int (x 11 , Xq) (A16) 

where the matrices Jj are the generators in the reducible representation of Spin(d) corre- 
sponding to ipint- The matrix elements of Ji are independent of x M : 

J d d x $1 (x m , x") 7, $nt (* m , xn = J d d x tilt (x"\ a#) Ji tint (x m , x£) . (A17) 

The primordial condensate is in a specific representation, but the basis functions in other 
representations are chosen to rotate with it according to ( 1A13I) and (1A15I) . 

It may be helpful to illustrate the above ideas by returning to the 2-dimensional analogy. 
Equation flA7j) becomes 



v (x) = 7l vec v (x Q ) , TZ vec = 
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or | 








\R(r) / 



for the vector representation and 

s(x)=n sp s(x ) ,n sp = e- l ^\ s(x )= \^~^ / J or ^ ~ I (A19) 

for the spinor representation. The matrices corresponding to the Jj are 

( -i \ a* 1 ( 1 \ 

J vec = and J sp =^ = - • (A20) 

yi o j 2 2 yo -l j 

Notice that fa is an angular coordinate in the internal space, whereas fa is a parameter 
specifying the rotation of iplnt a ^ fixed x m as x^ is varied. 



Appendix B: Example of solutions in internal space 



Our goal in this appendix is merely to show that there are solutions with the form required 
in Appendix |AJ so we will look for solutions with the higher-derivative terms in (I5.2p and 
(15. 3 p neglected, and with *f? int sufficiently small that Vq (x m ) can also be neglected. Then 
(15 .2p and f !5. 31) become 



2m 



dm^m i^int I ^ ' int (*£ j '-Ej 



2m 



C ( (^,/) = 0. (Bl) 



For simplicity of notation, let i\f int (x m , x^ 1 ) again be represented by i^{x), with com- 
ponents ip p (x). Each component varies with position in the way specified by ( 1A7h (to- 
gether with the radial dependence of ip(xo)). It therefore has a kinetic energy given by 
— (2m ) 1 dmdm'ipp (x), and an orbital angular momentum given by the usual orbital angu- 
lar momentum operators Jj in d dimensions 
ip p (x) varies as a function of the angles 0j. 

The Laplacian d m d m can be rewritten in terms of radial derivatives and the usual J 5 , 
giving 
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4461]. which essentially measure how rapidly 



42j-|44| 



1 d 



,2K 







J 2 



— J+--1 \^ p {x) 







K 



d 



(B2) 



after rescaling of the radial coordinate r. In addition, it is shown in Refs. 

J 2 ^ (x) =j(j + d- 2) i) p (x) 



Mkm that 



(B3) 
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where j is the orbital angular momentum quantum number, as defined on p. 677 of Ref. 43], 



but with this definition extended to half-integer values of m a and j. Normally, of course, 
only integer values of these orbital quantum numbers are permitted 4t|. However, the 
functions ip p (x) as defined in Appendix IA1 can have j = 1/2 etc. (in which case they are 
multivalued functions of the coordinates but single-valued functions on the group manifold, 
as discussed below (1A7j) ). Also, the demonstration of ( IB3j) in Ref. [43J can be extended in 
the present context to half- integer j, because it employs raising and lowering operators. (At 
each x, ip p is a linear combination of states with different values of m a , but ( 1B3[) still holds.) 
For each ip p (x) the radial wavefunction then satisfies 



1 d 



,2K 



d 

dr 



J U + d- 2) 



R (r) = . 



(B4) 



r 2K ^ r 

It may be helpful once again to consider the 2-dimensional analogy of Appendix |AJ where 
the orbital angular momentum operator is 



J = —id/d(f) . 



(B5) 



For the vector representation, (1A18I) implies that the kinetic energy is given by 



d m dmV (x) 



Id f d 



r dr \ dr 

\d_( d_ 
r dr V dr 



+ 



1 d 2 



*2 f)A\1 



cos <p — sin ■ 
sin <b cos d> 



V (£C 



v (x) 



in agreement with (1B4[) for j = 1. For the spinor representation, (1A19j) gives 



dm.dm.S (x) 



19/9 



r dr \ dr 



+ 



T*2 f)sf$ 



1 d f r d\ 1/4 



r dr V dr 



s (x 



in agreement with (1B4j) for j = 1/2. 



Equation ( 1B4I) can be further reduced to [45|, |46|] 

rf 2 k(k - 1 



dr* 



+ 



- 1 



X (r) = , k = j + K = j 



d-1 



(B6) 
(B7) 



(B8) 
(B9) 



(BIO) 



where x ( r ) = r R (r) . It is then easy to show that 



X (r) oc r k as r — y , x (r) oc sin (r + 5) as r — >■ oo 



(Bll) 
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a +-j-l)x(r)=0 (B12) 



where 5 is a phase. 

As a specific example, for d — 10 and j = 5/2 (i.e., the usual spinorial representation) 
we have 

\ dr 2 r 
with the exact solution 

X (r) = cX^sm^J (B13) 

21 630 18 900 496125 9 823 275 108 056 025 .„ x 
X = ! + ^ + — + ^^ + ^^ + ^T- + ^ ( B14 ) 

(where c is a normalization constant), which reduces to 

X (r) oc r 7 as r — )■ , x ( r ) ^ srn ( r + 5) as r — )• oo . (B15) 

Again, these examples are merely meant to indicate that solutions with the right basic 
form do exist. More realistic solutions will require detailed numerical calculations. 

Appendix C: Euclidean and Lorentzian Propagators 

For Weyl fermions, the Euclidean 2-point function is 



r ( \ 1 1 t mU A I V ^f V ^f^f ( X ^J M e ~ Sf rri x 
G, (x u x 2 ) = ( Xl )^ f (x 2 )) = Jv ^ v ^ f Js f ( C1 ) 

n s / # ; * (s) j # y (s) e-V M^W*) J2 suS2 ( Sl ) ^* ( S2 ) £/ (a*, s x ) C/t (x 2 , s 2 



(C2) 



n s / diP* f (s) J diP f (s) e -^f(sMWf(s) 

where (16 .2\\ and H6 .31) have been used. In a term with s 2 7^ S\, the numerator contains the 
factor 

/ chpf (si) / c% (s x ) e-^^i)^-!)^/^)^ ( Sl ) = (C3) 
according to the rules for Berezin integration. But a term with s 2 = si contributes 

j#;( 3 oj# / ( 3l )e-*;^.H^<-.»^( 3l )^( 3l ) t 

* (si) / (si) e^/W^O^O 

= a(s 1 )- 1 t/(x 1 , Sl ) ^(xz.si) (C4) 
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so 



Gf(x u x 2 ) = Y, a Gf(s)U(x u s)Ui(x 2 ,s) , G f (s) = a(s)- 1 . (C5) 

If the U (x, s) used to represent ipf (x) are a complete set, the propagator Gf (x, x') is a true 
Green's function: 

L f (x) U (x, s) = a (s) U (x, s) , ijj f (x) = J2 U s ) ~$f ( s ) ( C6 ) 

s 

and J2s U ( x -> s ) (x' : s) — S (x — x') imply that 

Lf(x)Gf (x, x') — 8{x — x') (C7) 

as usual. 

The treatment for scalar bosons is similar: 

g. (»„»») - w^fa)) = 7:f:f b f 2) e ~ s ' <«) 

\ / J V <j) ] b V <j) b e b f 

= UsS^d^MS^dlm^js) e -^)[(Hc? t W) a + (lm^W) a ] EsiiS2 fe ( Sl )0; ( 32 ) 

I!, H, d Re 6 (s) d Im 4> b (s) e~~ a{s) [(^(«)) a +(^«) a ] 

xf4 0r 1)Sl ) C/ 6 f (x 2 , S2 ) (C9) 

where 

L b (x)U b (x, s) = a(s)U b (x,s) , 6 (x) = ^ U b ( x , s ) 0& (s) ■ (CIO) 

s 

In a term with s 2 ^ Si, the numerator contains the factor 

dRe0 6 ( Sl ) [°° dlm^ b ( Sl ) e -5(^)[(R^ 6 ( S i)) 2 +(im^( Sl )) 2 ] [ Re b ( Sl ) + ,im0 b ( Sl )] 



= (Cll) 

since the integrand is odd. But a term with s 2 = si contains the factor 
/^dHe^CaOe- 3 ^^' 1 )) 2 (Re0 b ( Sl )) 2 /^dlm^ e -=(-0(im*(.i)) a (lm^( Sl )) 2 



= a(si) _1 (C12) 



so 



G 6 (m, x 2 ) = V G 6 (s) C/ 6 On, s) Ul (x 2 , s) , G b (s) = a (s)' 1 . (C13) 

< ' s 
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As usual, a (s) and a (s) contain a +ie which is associated with a convergence factor in 
the path integral (and which gives a well-defined inverse). 

The above are the propagators in the Euclidean formulation. The Lorentzian propagators 
are obtained through the same procedure with a (s) — > —ia (s) and a(s) — > —ia(s): 

C) (s) = ia (sy 1 , G L b (s) = ia (s)" 1 . (C14) 

The propagators in the Euclidean and Lorentzian formulations thus differ by only a factor 
of i. More generally, in the present picture, the action, fields, operators, classical equations 
of motion, quantum transition probabilities, propagation of particles, and meaning of time 
are the same in both formulations. (However, this equivalence holds only when there are no 
self-interactions of the form (^)\(j)b j , and ( 14. 5 p appears to imply that the original Euclidean 
formulation should be used when terms like these are significant.) 

For a single noninteracting bosonic field with a mass m&, the basis functions are 



U b (x, p) = V~ 1/2 e ip - x = y-i/2 e -^V^ (C15) 



so with s->pwe have 



and 



a (p) = to 2 — |~j7"| 2 — ml + ie (C16) 



G b (p) = - . 2 ^— (C17) 

lo z — \ p \ — ml + ie 

G\ (p) = - * ^— • (C18) 

u z — \ p \ — mi + ie 



Notice that ( 1C5[) and (1C13H hold even when the basis functions in ( 1C6[) or ( ICIOP are not 



a complete set. The propagators for the scalar bosons of Section IVlII therefore have the form 
of (1C13p . but include only those modes with s > 0, as specified in (I7.38p . This implies that 
radiative corrections will be modified. 
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